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ABSTRACT 

Aims. We introduce analytical response functions and their main properties as an important diagnostic tool that help understand Stokes 
profile formation physics and the meaning of well-known behaviors of standard inversion codes of the radiative transfer equation often 
used to measure solar magnetic fields. 

Methods. A Milne-Eddington model atmosphere is used as an example where response functions are analytical. A sample spectral 
line has been chosen to show the main qualitative properties. 

Results. We show that analytic response functions readily provide explanations for various well-known behaviors of spectral lines, 
such as the sensitivity of visible lines to weak magnetic fields or the trade-offs often detected in inversion codes between the Milne- 
Eddington thermodynamic parameters. We also show that response functions are helpful in selecting sample wavelengths optimized 
for specific parameter diagnostics. 

Key words. Radiative transfer, response functions, magnetic fields, spectropolarimetry, solar magnetism. 



1. Introduction 

Diagnosing the solar atmosphere from spectropolarimetric 
observations is one of the most important subjects of current 
solar physics. Both the theoretical understanding of the physical 
processes taking place in the photosphere and the design of 
new instrumentation that improve our ability to obtain more and 
better information from the Sun can be improved by a thorough 
study of the radiative transfer equation (hereafter referred to as 
RTE). RTE is the only tool available to describe the problem 
mathematically. Approximations have been devised according to 
the observational and the post-facto computational capabilities. 
The Milne-Eddington (M-E) approximation has provided insight 
into the physical processes taking place in line formation and 
inferring accurate values of several physical parameters of the 
solar atmosphere. Its specific analytical character is its most 
powerful feature. 

A physical analysis of the sensitivities of spectral lines in 
terms of analytic mathematical functions is still missing in the 
literature and may provide a better understanding of how the 
solar parameters influence the shape of the observed Stokes 
profiles of these spectral lines and explanations for the trade- 
offs and other well known behaviors of inversion codes currently 
used for the inference of such solar atmospheric parameters. 
Here we introduce the analytic response functions (RFs) of 
Stokes profiles as formed in M-E model atmospheres and discuss 
their main properties. 

Weighting functions for unpolarized light (Mein 1971) were 
the precursors of RFs, extended to polarized light by Landi 
DeglTnnocenti & Landi DeglTnnocenti (1977). As explained 
by Ruiz Cobo & del Toro Iniesta (1994; see also del Toro 
Iniesta 2003), RFs provide the sensitivities of Stokes profiles 
to the various atmospheric quantities playing a role in line 
formation. Since all these quantities are constant with depth 



in a M-E atmosphere, M-E RFs are simply partial derivatives 
of the analytic solution of the RTE with respect to the model 
parameters. This feature enables us to deduce analytic formulae 
for the sensitivities (they are explicitly written in the Appendix) 
and to study their characteristics and properties. Such properties 
are shown to be useful in practice in understanding the behavior 
of spectral lines as well as in helping in line and sample selection 
when designing new instruments. 

2. The response functions in a Milne Eddington 
atmosphere 

2.1. Summary of radiative transfer 

The radiative transfer equation (RTE) for polarized light in a 
plane-parallel atmosphere is 



f=K(I-S), 
dr 



(1) 



where / = (/, Q, U, V) T stands for the Stokes vector which gives 
a full description of the polarization state of light, t for the 
continuum optical depth at a reference wavelength, K for the 
4x4 propagation matrix, S for the source function vector and 
T means the transpose. All the medium properties relevant to 
line formation are contained in K and S. In LTE conditions, 
S = (B A (T), 0, 0, 0) T , where B A (T) is the Planck function at the 
local temperature T. 

In a Milne-Eddington (M-E) model atmosph ere, an ana- 



lytical soluti o n is f o und fo r the RTE (see, e.g. Unnol 19561: 



Rachkovskvl I1962L 1 19671 iLandolfi & Landi DeglTnnocentil 
19821) . In such an atmosphere, all the atmospheric quantities are 
constant with depth except for the source function that varies 
linearly: 
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(2) 
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The propagation matrix is also constant wit h depth. 
Following, e.g., the notation in Idel Toro Iniestal (120031) . such an 
analytical solution reads 

/ = S o + A -1 [77/(77? + p 2 Q +pl+pl)] S 1, 

Q = -&~ l bfiT]Q + TiiOivPu -vuPv) + Pell] Si, 

U = -A^tftfu + tu(t]qp v - T] V p Q ) + puYl]S u 

V = -A -1 [rfjriv + t]i(VuPq ~ VqPu) + Pv^] Si, (3) 

with 

A = rjjirfj -rf-Q-rjl-rfy+pl+pl + p\) - IT 2 , (4) 
where 



II = 77 2 p e + TjuPu + TJvpy. 



(5) 



It can easily be seen that 77/, tjq, r]u, rjy, pq, pu, and pv, 
and hence the solution depends on just nine parameters, (B, 
y, x)> the strength, inclination and azimuth of the magnetic 
field vector on the local reference frame, on So, Si, the two 
parameters describing the source function, on 770, the line-to- 
continuum absorption coefficient ratio, on A/Id, the Doppler 
width of the line, on the damping parameter a, and on the line- 
of-sight velocity, vlos- 

2.2. Milne-Eddington response functions 



According to iRuiz Cobo & del Toro Iniesta (1 19941) (see also 



Idel Toro Iniesta & Ruiz Cobol 1 19961 or Idel Toro Iniestal 120031) 
the sensitivity of the Stokes profiles to perturbations of the 
atmospheric physical quantities is given by the response 
functions (RFs). Fortunately, in the specific case of constant 
quantities (model parameters) with depth, as is the case of an 
M-E atmosphere, such RFs are the result of integrating in depth 
the regular RFs. Such T-integrated response functions are thus 
simply functions of the wavelength and can be considered as 
the partial derivatives of the Stokes vector with respect to the 
corresponding model parameter: 



dig) 

dx 



(6) 



where x represents any of the model parameters. We hereafter 
refer to these r-integrated RFs as M-E RFs or just RFs. 

Therefore, by taking derivatives of the analytical solution 
(01, the sensitivities of the Stokes profiles to perturbations of 
the M-E model parameters can be found (see the Appendix for 
explicit formulae). These sensitivities are the only tools we have 
to evaluate our ability of determining the various quantities: 
should the / Stokes vector not vary after a perturbation of 
a parameter, x, we would be unable to infer it from the 
observations (it would not be a proper model parameter). 

2.3. Line sensitivities: the shape of M-E RFs 

Equations (01 and © provide the necessary means for studying 
the behavior of the M-E Stokes profiles. The shapes of the 
RFs do not vary dramatically from model to model or from 
line to line. M-E RFs appear homologous to each other. This 
property allows us to choose a single line to illustrate the 
practical usefulness of our functions. Let us take the Fe 1 line 
at 525.064 nm as an example. We sele ct this line because 
it is used by the IMaX magnetograph ( iMartmez Pillet etail 
12004 and some of the results have implications either for the 



design or for the analysis of the data to be obtained with this 
magnetograph. The line has an effective Lande factor of 1.5 
and is often consi dered to be quite in sensitive to temperature 
perturbations (e.g.. IStenflo et al.l [l984). A single model is also 
enough for our purposes. We have used the NSO Fourier 
Transform Spectrometer atlas as a reference spectrum and the 
line was fit with errors smaller than a 2%. The resulting model 
parameters are: So = 0.02, Si = 1 , 770 = 7.2, a = 0.3, A/Id = 30 
mA and a macroturbulent velocity, v mac = 0.37 km/s. Unless 
otherwise stated, all the numerical examples that follow refer to 
this line and this model. Several magnetic field strengths (200, 
800, 1400 and 2000 G) have been used to synthesize the Stokes 
profiles and their RFs, assuming a constant field inclination and 
azimuth of 45 degrees. 

Fig.Q~]shows the synthesized Stokes profiles. As the magnetic 
field increases, the Stokes V lobes increase but their peaks 
do not separate much because the strong field regime has not 
been reached for this line with these strengths. In Fig. [2] we 
give a graphical illustration of the analytical RFs of the four 
Stokes parameters to magnetic field strength perturbations. Both 
the Stokes profiles and the RFs present wavelength symmetry 
properties, as expected from a M-E model atmosphere. The 
RFs to the magnetic field strength preserve the Stokes profile 
symmetries while velocity RFs display opposite parity (see 
Fig.EJ. 

Fig. 12 shows that the response of the line is wavelength de- 
pendent. Different wavelength positions have different sensitivi- 
ties. Within a single Stokes profile different wavelength samples 
react differently to the same perturbation. Some of the samples 
are insensitive. For instance, in this example the Stokes V zero- 
crossing point remains the same regardless of B and, hence, the 
response is zero at this wavelength. All the RFs show peaks cor- 
responding to different maxima and minima. Note that these ex- 
trema pinpoint where the Stokes profiles are more sensitive to 
perturbations of the physical quantity: the greater the peak, the 
greater the sensitivity. 

Although Stokes /, Q and U are more sensitive to B 
perturbations when the strength is greater, the Stokes V profile 
sensitivity to field strength perturbations is a maximum for the 
weak fields and decreases with increasing field strength. In the 
weak field regime, Stokes V is proportional to B and any change 
of B is translated directly into an increase (or a decrease) of the V 
signal; when the field increases, however, a competition between 
increasing the profile and peak separation becomes important; 
At a given B value, peaks will no longer increase but separate 
from each other. This behavior of Stokes profiles is known for 
long but a glance to the Stokes V panel of Fig. 2 illustrates it in a 
very clear way. Moreover, the sensitivity of Stokes V in the weak 
field regime explains the reasonably accurate inversion results 
for weak magnetic fields obta ined in numerical experiments by 
IWestendorp Plaza et al.l (Il998l) . 

Fig.|3]shows the Stokes RFs to LOS velocity. The first clear 
feature in this figure is that neither the sizes nor the shapes 
depend on the LOS velocity. The latter only shifts the RFs as 
it does with the profiles. The RF size is larger for Stokes / and 
V than for Stokes Q and U, because of the corresponding size of 
the profiles. Since Stokes / and V are larger than Stokes Q and 
U in this example, most information on velocities is carried by / 
and V. The LOS velocity can always be well determined because 
the loss of sensitivity to vlos perturbations of the Stokes / profile 
is compensated by that of the V profile. 

The Stokes / RF to LOS velocity perturbations decreases 
with B when the Stokes Q, U, and V RFs increase. This is due 
to the different shape ratios of the various profiles. According 
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Fig. 1. Stokes I/I c , Q/I c , U/I c and V/I c , for the Fei line at 525.06 nm, with a magnetic inclination and azimuth of 45 degrees. 
Different lines stand for different magnetic field strength values. The Stokes parameters are normalized to the local continuum. 



to Cabrera Solana, Bellot Rubio and del Toro Iniesta (2005), the 
spectral line sensitivity to the LOS velocity is mostly determined 
by the ratio between the width and the depth of the line. The 
greater the field strength, the wider and shallower the Stokes / 
profile. Therefore, its sensitivity to vlos perturbations decreases 
with increasing B. Each lobe of Stokes V, however, first becomes 
larger and then narrower and steeper at the central wavelength as 
B increases. Hence its greater sensitivity to vlos for increasing 
field strengths. 

The relative maxima of the RFs to LOS velocity perturba- 
tions correspond to wavelength positions where the inflection 
points of the Stokes profiles are located independently of the 
model atmosphere and spectral line. For instance, the minimum 
of Stokes / and the peaks of Stokes V correspond to zeros on the 
corresponding RFs to LOS velocity, therefore regions where the 
Stokes profiles do not change when LOS velocity does. 

The extrema of the RFs to B and to vlos perturbations do not 
coincide with those of the corresponding profiles. This fact can 
be clearly seen in, e.g., the bottom right panels of Figs. [2] and 
[3] Therefore, the extrema of the Stokes profiles do not carry, 
in principle, more information on given parameters than any 
other particular wavelength sample. Another very interesting 
feature is that, for a given spectral line, the RFs differ from 
each other. RFs to magnetic field strength perturbations do not 
resemble those to LOS velocity perturbations (compare Figs. [2] 
and [3]). For instance, their maximum sensitivities (RF peaks) 
are placed at different wavelengths. These differences among 
RFs help disentangle the influences on spectral line formation 
of the various model quantities and allow inversion algorithms 



based on RFs (e.g., SIR by iRuiz Cobo & del Toro Iniestal 
Il992h to obtain accurate results: if a given Stokes profile is 
inappropriate for a particular wavelength sample, other profile 
or wavelength samples provide the required information. RF 
differences can also be seen for the other M-E parameters 
except for A/Id, ??o an d a - The RFs to these thermodynamic 
parameter perturbations are very similar to each other as can 
be seen in Fig. [4] A small perturbation of any of these three 
parameters produces a modification in the Stokes profiles that 
is very similar to the changes produced by small perturbations 
of the other two. These similarities between the A/Id, 77o and 
a RFs explain the trade-offs often observed in M-E inversions. 
Fortunately, their RFs are different enough from those of the 
other mo del parameters for them to be accurately retrieved 
(see, e.g. JWestendorp Plaza et alj|1998l) . Thus, the M-E model 
atmosphere, although providing a simplistic scenario for line 
formation which may not full account for all thermodynamic 
properties, allows fairly accurate inferences of the constant 
magnetic field vector B, vlos, Sq and S i . 

The RFs to magnetic field inclination and azimuth 
perturbations do not depend on the derivatives of the absorption 
and dispersion profiles; thus, the shapes of the RFs are very 
similar to the corresponding Stokes profiles (see Fig. [5]). Only 
Stokes Q and U respond to azimuth perturbations. The larger the 
field strength, the greater the sensitivity of the Stokes profiles to 
y and x perturbations. This is again an explanation of a well 
known fact: we measure y and x more accurately when B is 
strong. 
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Fig. 2. Analytical M-E RFs of Stokes I/I c , Q/I c , U/I c and V/I c to magnetic field strength for the Fei line at 525.06 nm, with a 
magnetic inclination and azimuth of 45 degrees. Different lines stand for different magnetic field strength values. Units are 1CT 3 
G '. 



2.4. Relative response functions 

So far we have only discussed "absolute" RFs, i.e., functions 
with dimensions; e.g. the RF to B is measured in G , that 
to vlos perturbation is measured in (km s -1 ) -1 and so on: 
RFs give modifications of the profile per unit perturbation 
of the parameter. To co mpare them to one another, r elative 
RFs should be used dRuiz Cobo & del Toro Iniestal Il994t 



Idel Toro Iniesta & Ruiz Cobol 119961) . These relative responses 
are obtained by multiplying the standard, absolute RFs by 
the corresponding model parameter. Relative RFs tell us how 
sensitive one model parameter is compared to the others. For 
instance, the relative RF to A/Id is much larger than that to f/o and 
that to a (in particular three times as large as the RFs to 770 and 
twenty times larger than those to a for Stokes /, in our sample 
M-E atmosphere). This means that a small relative perturbation 
of A/Id changes the Stokes profiles much more than the same 
relative perturbation of ?/o or a. Consequently, A/Id should be 
better determined by M-E inversion codes. 

2.5. Two-component model atmospheres 

Model atmospheres with two or more components are 
commonly used in the analysis of observations. Any two- 
component model atmosphere is based on the assumption that 
within the resolution element two different atmospheres coexist, 
namely, one magnetic atmosphere filling a surface fraction a, 
and one non-magnetic atmosphere filling the remaining (I - a) 



fraction, a is called the magnetic filling factor. If I m stands for 
the Stokes profile vector emerging from the magnetic region and 
/ m for that of the non-magnetized atmosphere, the observed 
Stokes vector can be written as / = I m a + 7 nm (l - a). 

Thus, according to Eq. (O, the RFs to a perturbations are 
given by I m - I nm . Hence, the larger the difference between 
the magnetic and the non-magnetic atmospheres, the greater the 
sensitivity to a. Since most of the difference is the polarization 
signal itself, Q m , U m , V m , when this signal is strong we can 
easily discern it from the non-magnetic signal. 

2.6. The influence of smearing 

Spectral line smearing by macroturbulence is a well known 
effect that needs be taken into account in the analysis of most 
observation s except, perhaps, in those with very high spatial 
resolution dAsplund et al.l l2000h . Besides macroturbulence, 
instruments have finite-width profiles that produce smearing 
of the observed Stokes profiles which become wider and 
with smaller peaks. This smearing reduces the information 
on physical parameters carried by the spectral line through 
convolution: / t, s = / * F(A), where * stands for the convolution 
symbol and the scalar smearing profile, F(A), is convolved with 
all the four Stokes parameters. 

This loss of information through smearing is also translated 
into a loss of sensitivity of spectral lines to the atmospheric 
quantities. In fact, since the derivative of a convolution is equal 



4 



D. Orozco Suarez and J.C. del Toro Iniesta: The usefulness of analytic response functions 



> 
o 

-4 — 1 

to 
o 

-t—> 

CO 
Li. 

cr 




-200 -100 m 100 
Lambda (mA) 



200 



> 
o 

-t—s 

a 

in 
<u 

o 

-t—s 

in 

U. 

cr 



0.004 
0.002 
0.000 
-0.002 
-0.004 



/ \ 




-200 -100 J 00 
Lambda (mA) 



200 



g 0.04 



> 
o 

ZD 

in 
a> 

o 

on 

Li_ 

OH 



0.02 
0.00 
-0.02 
-0.04 



. 1 ■ ■ 


1 

/\ 


i i . 




/' \ 






















— — — 






>,*••.» 






\V 


\\ I- 






»w.' 






\ 1 


■ 1 . . 


■ 





-200 -100 < 100 
Lambda (mA) 



200 



> 
o 



in 
o 



0.10 
0.05 
0.00 



w -0.05 

u. 

oh 

-0.10 



J 1 

- //.•'•• 


1 1 1 

200 

^\ 800 

,\ 1400 . 

\'( 2000- 


/'•/\ 
/'/ 

■ ^ — ■// 


\\\ 


^> v\ V 
\ \ •• V 

V 


If/// 


' 





-200 -100 



100 



200 



Lambda (mA) 



Fig. 3. Analytical M-E RFs of Stokes I/I c , Q/I c , U/I c and V// c to LOS velocity for the Fei line at 525.06 nm, with a magnetic 
inclination and azimuth of 45 degrees. Different lines stand for different magnetic field strength values. Units are [km/s] ■ 



to the convolution of the derivative of one of the functions with 
the second one, response functions become smeared as well: 

Robs,, = R x * F(A). (7) 

Fig. [6] shows the effect of RF smearing. The convolved RFs 
are smoother and significant information is lost. 

2. 7. The influence of noise 

Stokes profiles are affected by the noise intrinsic to the 
observational process. Should the polarimetric signal be buried 
by noise, any algorithm one could devise to determine 
atmospheric quantities would fail. Therefore, our ability to infer 
accurate solar parameters depend significantly on the signal-to- 
noise ratio of the observations. 

Response functions can help in quantifying this effect. Since 
RFs simply provide the modification of the Stokes profiles after 
a perturbation of the physical quantities, if that modification is 
smaller than the noise level it will be effectively undetectable. In 
other words, the size of RFs to perturbations of a given quantity 
sets a threshold for the detection of a unit of such a quantity: 
for instance, according to Fig. [21 1 Gauss will only be detectable 
by a single wavelength sample if the noise is below 1.5 TO -4 
(continuum intensity is at 1); within the linear approximatiorfl, 
10 Gauss will be detectable with a noise below 1 .5- 10 3 and so 
on. 

1 RFs come from a linear perturbation analysis of the radiative 
transfer equation 



According to Eq.( [6]), a standard deviation in the Stokes 
signal cr will induce an error cr x per single wavelength sample 
given by: 

* = 4 R lx +R l* + R lx + R l* a *- (8) 

Detectability should increase, of course, as the root of the 
wavelength sample number, but the RFs allow an estimate of 
the expected accuracy of our inferences. 

The above estimates can be considered lower limits for the 
errors since model parameters are not independent of each other 
and correlations may exist between sensitivities such as those 
already reported between 770, A/Id and a. 

3. The usefulness of the RFs for vector 
magnetographs 

Modern vector magnetographs are not restricted to using one 
or two wavelength samples as are the classical magnetographs. 
Instruments like IMaX are devised to measure up to five 
wavelengths: one in the continuum and four across the line 
profile. The choices of the spectral line, of the number of 
samples and of the precise wavelength for each of them are 
important issues that arise during the design phase of the 
instrument. This section is aimed at illustrating how the RFs can 
help such decisions. 

Finding a suitable spectral line is crucial and can be 
ach ieved through RFs on the sim ple phenomenological model 
by ICabrera Solana et"aT1 d2005l) that allows establishing a 
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ranking of sensitivities to the different atmospheric parameters 
among the various lines considered. The IMaX Fei line at 
525.064 nm can be seen in Fig. [8] Da ta for this line have 
been i ncluded in the original figure by ICabrera Solana et"aT] 
(120051) . where it is identified as one of the most sensitive 
of the set to velocity perturbations. It has a medium 
sensitivity to magnetic field strength perturbations in both 
the strong and the weak field regimes. However, it is not 
very sensitive to temperature (not shown) and, hence, a 
good candidate for inferences in the various solar structures 
avoiding thermodynam ic trade-offs. The Helioseismic and 
Magnetic Imager (HML lScherrer & SDO/HMI Teamll2002l) and 
the V isible-light Imager and Magnetograph (VIM: lMarsch et"a"fl 
120051) . two planned instruments for the Solar Dynamics 
Observatory, NASA, and the Solar Orbiter, ESA, missions, will 
use the Fei line at 617.334 nm. This spectral line is very well 
ranked in Fig. [8] for inference of both magnetic field strengths 
and LOS velocities. 

A minimum number of wavelength samples is obtained by 
roughly doubling the free parameters of the model: since an M- 
E model is made up of just ten parameters, a minimum of twenty 
observables (five wavelengths times the four Stokes parameters) 
is needed. This is the choice for all the three instruments 
mentioned above. Unfortunately, no purely objective means 
exists to select the wavelengths for the samples. Nevertheless, 
RFs are a powerful tool that help select those wavelengths that 
better suit our purposes. If one is interested, for instance, in 



just the magnetic field strength and neglects the other physical 
quantities, choosing those wavelengths where the RFs to B 
reach local maxima would be appropriate. If the interest lies in 
several physical quantities at the same time (e.g. on the three 
components of the magnetic field and on the LOS velocity) we 
suggest the use of a linear combination of regular RFs weighted 
according to the specific interests. Since RFs can be positive 
or negative, we propose the use of absolute-valued RFs. Hence, 
consider 

(9) 

i 

where j runs from 1 through 4, corresponding to the four Stokes 
parameters, and index i accounts for the physical parameters. 
Since the set of weights a, can be tailored at will, there is no 
single choice for samples but an examination of ft provides 
important hints for the selection. As an example, Fig. Q shows 
different linear combinations for the IMaX line case. If index i 
runs from 1 through 4 standing for B, y, (p, and vlos, respectively, 
the plotted curves correspond to ££1,2,3,4 = 1,1,1,1 (solid, 
black lines), 121,2,3,4 = 2,2,2,0.5 (dashed, blue lines), 0:1,2,3,4 = 
3, 1, 1,0.5 ( dashed-dotted, red lines), and ai,2,3,4 = 0.5, 1,1,3 
(dotted, green lines). The vertical lines indicate a possible choice 
for wavelength sampling (+42, 78 mA), selected mostly from the 
properties of the Stokes / and V RFs since these two parameters 
usually exhibit the largest signals in solar atmospheres. While 
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Fig. 5. Analytical M-E RFs of Stokes I/I c , Q/I c , U/I c and V/I c to magnetic field inclination, y, for the Fe i line at 525.06 nm, with a 
magnetic inclination and azimuth of 45 degrees. Different lines stand for different magnetic field strength values. Units are in 1CT 3 
[degrees] -1 . 
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Fig. 6. Analytical M-E RFs of Stokes V/I c to the magnetic field strength (left panel) and of Stokes I/I c to the LOS velocity (right 
panel) for the Fe i line at 525.06 nm, with a magnetic field strength of 2000 G and field inclination and azimuth of 45 degrees. The 
dashed lines correspond to the RFs convolved with a Gaussian smearing profile of 60 mA of full width at half maximum (FWHM). 
Solid lines correspond to the original RFs. Units are in 10~ 3 G~' (left) and [km/s] -1 (right). 



the most external samples seem to be optimum, some other good 
choices for the inner wavelengths are possible and up to the user. 



4. Conclusions 

The many interesting features of analytic response functions 
have been discussed in this paper by considering the specific 
case of an M-E model atmosphere. Since an analytic solution for 
the radiative transfer equation is available for this atmosphere, 
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Fig. 7. Four different linear combinations of the Stokes vector RFs for the IMaX line. The plotted curves correspond to 
<*i,2,3,4 = 1, 1, 1, 1 (solid, black lines), 0:1,2,3,4 = 2,2,2,0.5 (dashed, blue lines), «i,2,3,4 = 3, 1, 1,0.5 (dashed-dotted, red lines), 
and ai,2,3,4 = 0.5, 1, 1,3 (dotted, green lines). The light-grey, vertical lines indicate a possible choice for wavelength sampling 
(±42, 78 mA). 



the sensitivities of spectral lines, as described by RFs, can also 
be cast in an analytical form by simply taking partial derivatives 
of such a solution with respect to the model parameters. The 
analytic M-E solution has been thoroughly used in the past for 
insight into radiative transfer physics and as a purely practical 
diagnostic tool through the M-E inversion codes. Likewise, we 
have shown in this paper that the analytic, M-E RFs are useful 
to better understand spectral line formation and the behavior of 
Stokes profiles in different formation conditions and also for 
practical recipes that can help in selecting spectral lines for given 
purposes, in selecting wavelength samples, etc. 

A summary of the various results obtained follows: 

1. Response functions look homologous to each other, hence 
enabling qualitative, general discussions by considering a 
single spectral line in a specific model atmosphere. Here, we 
have targeted the IMaX, Fe 1 line at 525.60 nm in a M-E 
model representative of the quiet Sun thermodynamics (as 
observed by FTS) and with various vector magnetic fields 
and LOS velocities. 

2. The sensitivities of spectral lines to the various parameters 
depend on wavelength: some samples are better suited to 
diagnose a given parameter; some wavelengths are even 
insensitive to another parameter. The RF extrema show 
trivially those wavelengths where sensitivity is maximum. 

3. As expected in M-E conditions where no gradient of LOS 
velocity is present, RFs display clear wavelength symmetry 



properties. The RFs to magnetic field strength perturbations 
show similar parity as the Stokes profiles while the RFs to 
LOS velocity perturbations are of opposite parity. 

4. Stokes V sensitivities to B perturbations are significant for 
very weak field strengths. This fact explains the reasonably 
accurate results of M-E inversions in this strength regime. 

5. The shape of the RFs to LOS velocity perturbations does not 
depend on vlos except for the Doppler shift. Variations of 
sensitivity of the Stokes / and V profiles are compensated: 
when information on vlos decreases in Stokes / it increases 
in Stokes V, so that vlos remains well inferred in any 
circumstance. 

6. We understand the trade-offs often found in the inversion 
codes among M-E thermodynamic parameters: their corre- 
sponding RFs are very similar to each other. Fortunately, 
they are different from the other RFs and can accurately in- 
fer vector magnetic fields and LOS velocities. Among the 
thermodynamic parameters, the relative sensitivity to AAd 
perturbations is larger than that to rjo and a, hence enabling 
better inferences. 

7. The magnetic filling factor a is better determined if there are 
significant differences between magnetic and non-magnetic 
atmospheres. When B is large this result is natural; when 
B is small, this result explains that differences in the 
thermodynamics of both atmospheres can help in inferring 
a properly. 
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8. Direct estimates of affordable noise levels can be directly 
obtained from RFs. 

9. Response functions can also be used to select spectral lines 
for given purposes or for given measurements. Moreover, a 
suitable combination of RFs provides quantitive arguments 
for wavelength sample choice. 
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Appendix A: Explicit formulae 



The propagation ma trix K of the RTE can be cast in the form (e.g. 
del Toro Iniesta 2003) 



K : 
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(A.2) 



and <t> p ,b,r and ifr p j v are the absorption and dispersion profiles, the p, b, r indices 
stand for the n and <r components of a Zeeman multiplet, and i]q is the ratio of 
line to continuum absorption coefficients. 

<p Pt b,r and tppju can be written as a sum of as many absorption and 
dispersion profiles as the number of p, b, r components as follows: 

<t>j = S M ,M„.jH(a,v), 

M,-M„=j 

= 2 ^ S M,M,„jF(a,v), (A.3) 

M,-M„=j 

S M[M„ j being the strength of each component with j = -1,0,1 corresponding 
to b, p and r. v stands for the wavelength in Doppler units which follows 
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The wavelength shift of the different Zeeman components with respect to 
the original position is given by 



A,l B 



eAlB 
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(A.5) 



4nmc 2 

where / and u stand for the lower and upper levels of the line transition, g for the 
level Lande factor, and M for the magnetic level quantum number. 

The evaluation of RFs reduces to the derivatives of the Stokes 
vector, / = (/, Q, U, V), with respect to the nine parameters, 
(Bo, B\, rjo, B,y,x, A/Id, ^los, o). In order to easily show such derivatives 
suppose a generic parameter x. Then, 
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where for simplicity 
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A and n are defined in Eqs. (4) and (5), respectively. Their derivatives are 
thus given by 
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The derivatives with respect to m can be easily calculated from Eq. (A. 2.) 
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The derivatives with respect to the other parameters imply the derivatives of 
the absorption and dispersion profiles and these lead us to obtain the derivatives 
of the Voigt and Voigt-Faraday functions (as defined by Landi degl'Innocenti, 
1976): 
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By using the chain rule and the derivatives of H(a, v) and F(a, v) with to 
respect a and v, 
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Fig. 8. Upper panel: Maximum value of the r-integrated RF 
to vlos for different lines as a function of the shape ratio 
multiplied by the central wavelength of the transition (see 
ICabrera Solana et"aD 120051) . Middle panel: Maximum values 
of the integrated RF to B for the same set of lines with 
geff^O, as a function of the shape ratio multiplied by the 
squared central wavelength (strong field regime). Bottom panel: 
Maximum values of the integrated RF to B (weak field 
regime). The sensitivities have been evaluated in the quiet Sun 
(crosses), penumbral (circles) and hot umbral (filled circles) 
model atmospheres. Dotted, dashed, and dash-dotted lines mark 
specific transitions in the quiet sun, penumbral, and umbral 
models, respectively. 
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